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ABSTRACT 



Context. Previous authors have reported the detection of intrinsically faint sub-mm emission lensed by the cluster MS045 1.6-0305. They 
suggest that this emission arises from a merging system composed of a Ly-break galaxy and a pair of extremely red objects which are 
multiply-imaged in the optical/NIR observations. 

Aims. Since the submm emission presents an unusually large angular extent (~ 1'), the possible radio emission asociatted with that system can 
help to identify optical/NIR counterparts due to the higher spatial resolution and astrometric accuracy of the radio observations. 
Methods. Archive VLA data (BnA configuration at 1.4 GHz) was reduced and analysed. A simple lens model was constructed to aid the 
interpretation of the radio and pre-existing sub-mm and optical/NIR data. 

Results. We present a 1.4 GHz map of the central region of MS045 1.6-0305 and report the detection of gravitationally lensed radio 
emission, coincident with the previously discovered sub-mm lensed emission. The overall morphology and scale of the radio and sub-mm 
emission are strikingly similar, extending ~ 1' across the sky. This observation strongly suggests that the radio and sub-mm emission 
arise from the same sources. Preliminary estimates of the total Sgso/jm/S iaghz flux density ratio appear to be consistent with that expected 
from distant star forming galaxies. The radio emission is resolved into 7 distinct components, and the overall structure can be explained, 
using a simple lens model, with three multiply-imaged radio sources at z ~ 2.9. One of these sources is predicted to lie in the middle of the 
previously mentioned system in the source plane, suggesting that it is related to the intense star formation generated during the merging process. 
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1. Introduction 



Sub-mm galaxies (SMGs) were first detected by scubaQ 



jSmail et alJI 19971) and are beli eved to be dusty star forming 



galax i es located at high reds hift jlvison et al .l2002t l ISmail et al 



2002t IChapman et alj|2003h . It is also suggested that they are 



the progenitors of present-day massive elliptical galaxies (e.g. 
Lillv et alj|l999l: ISwinbank et al.ll2006l) . Little is known about 
the objects associated with the faint end of the SMG popu- 
lation (SiiiO/jm < 2mJy), but they are pre dicted to domin ate 
(energetically) the population as a whole dKnu dsen 2004). A 
recent statistical stacking analysis ( iKnudsen et al.ii2005i) . sug- 
gests that distant red galaxies (DRGs) and Extremely Red 
Objects (EROs) contribute ~ 50% of the flux density of sub- 
mm sources with 0.5 < Sg,50fim < 5mJy. 

Intrinsically faint SMG cannot easily be detected, since 
their flux densities lie below the ~ 2 mJy confusion limit of 
SCUBA images at 850 fim. Typically they also fall well below 
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' Submillimetre Common-User Bolometer Array, mounted at the 
James Clerk Maxwell Telescope (JCMT) 



the (thermally hmited) sensitivity of current radio instruments, 
such as the VLA. Individual systems can often only be detected 
via strong gravitational lensing effects, produced by massive 



foreground clusters of g alaxies jKnudsenl I2004t iKneib et al 
20041 iGarrett et al.ll2005l) . 



The spatial magnification provided by the lensing cluster 
overcomes instrumental confusion limitations in the sub-mm 
and also boosts the measured flux density of the source (pro- 
vided the lensed images remain unresolved) , thereby increas- 
ing the probability of detection jBlainl Il997h . Another advan- 
tage of cluster lensing is that the magnification provided by the 
lens effectively increases the spatial resolution of the observa- 
tions, with the largest magnifications usually occurring in cases 
of multiple imaging. 

SMM J16359-I-6612, associated with the cluster Abell 
2218, was the first intrinsically faint, multiply i maged 



SMG dete cted in both the sub-mm ( iKneib et al.1 120041) 
and radio dGarrett et al.l l2005h . In this paper, we present 
VLA 1 .4 GHz radio observations of a second case. 



SMM J04542-0301 (IChapman et all l2002al) . associated with 
the cluster MS045 1 .6-0305. 
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Fig. 1. The VLA 1.4 GHz naturally weighted contour map (solid white lines) superimposed upon the SCUB A 850-//m contou r 
map (solid thin black lines) and the inverted HST F702W image of the centre of the cluster MS045 1 .6-0305 dBorvs et alj|2004 . 
The axes represent the right ascension (x-axis) and declination (y-axis) in the J2000 coordinate system. The solid thick black 
cu rves are the tangent ial (outer) and radial (inner) critical lines at z=2.9 1 1 associated with the lens model of the cluster determined 
by lBorys et al. The boxes are the positions obtained via Gaussian fits of the radio sources. The diamonds are the positions 

of three EROs from lTakata et al.l ( 2003 ). and the crosses/circles are the positions of a LBG lensed as tw o arc (ARC! and A RCl 
ci), and a triply-imaged EROs pair: B1/B2/B3 (images of ERO B) and C1/C2/C3 (images of ERO C) (see lBorys et"aDl2004 . The 
squares and circles have a size of 1" to illustrate the random and systematic errors due to measurement indeterminations and the 
aligment of the different images. Contours of the radio map are drawn at -3, 3, 4, 5, 8, 12 & 16 times thel-cr noise level of 9 /zJy 
per beam. Contours of the sub-mm map are drawn at 4, 6, 7, 9, 10, 11 & 11. 5 mJy per beam. The white circle inside a box in the 
bottom-right corner is the beam-size of the radio map (6.99 x 6.03 arcsec in position angle PA - 32.6°) whereas the black one 
corresponds to the beam-size of the sub-mm map (15 x 15 arcsec). 



MS045 1 .6-0305 is a cluster of galaxies situated at z=0.55 
( Gioia & Luppind 1994 ) that has rec ently been studied using 
optical and near-Infrared (NIR) data (iBorys et al.ll2004l) . They 
conclude that the sub-mm emission is probably related to an 
interacting system of three objects lying at z ~ 2.9: a Lyman 
Break Galaxy (LBG) and a pair of Extremely Red Objects 
(EROs). In the optical and NIR images, it is proposed that the 
LBG is imaged into two visible arcs, and the ERO pair are 
responsible for 5 additional sources of emission in the field 



dBorvs et alj|2004l : iTakata et alj|2003l) . However, the emission 



coming from the north-eastern and the central regions of the 
sub-mm image, cannot be fully rep roduced using the LBG and 
the ERO pair alone (see Fig. 7 from lBorys et al ] |2004 . 

In this paper, we present deep, high resolution 1.4 GHz 
VLA observations of SMM J04542-0301. In Sect. [21 we de- 
scribe the VLA data analysis and present the associated radio 
images. Section[3]describes a simple lens model for the system, 
in an attempt to explain the lensed nature of the radio emission 
related to SMM J04542-0301. In Sect.m we compare the ra- 
dio and sub-mm emission, including a discussion about possi- 
ble optical/NIR counterparts, and the preliminary calculation of 
the 5 sso^^m/^ 1.4GHz Aux density ratio. A summary of our main 
results is presented in Sect. |5] In the following discussion, we 
assume a ACDM cosmological model with Q.,„-Q3 , Q.x-Q.1 
and /jo=0.7 



2. Radio Observations 

VLA 1.4 GHz observations of the cluster MS045 1.6-0305 
were made in June 2002, and were retrieved from the NRAO 
data archive systerrQ. The integration time was 7.8 hours with 
the VLA in BnA configuration, employing two 25 MHz IFs 
in both left and right-hand circular polarization. Each IF was 
subdivided into 7 channels. The data analysis was performed 
using the NRAO AIPS package using standard analysis tech- 
niques. The absolute flux density scale was set by observations 
of 0137+331, and phase calibration was performed via short 
observations of 0503+020 between the 1 hour target scans. A 
wide-field image was made and bright sources far from the field 
centre were subtracted from the data. Self-calibration using the 
remaining sources in the centre of the field realised images with 
a Icr rm.s. noise level of 9 //Jy/beam. 

In Fig[T] we present the radio contour map (solid white 
lines) of the naturally weighted VLA image of SMM 
J04542-0301, superimposed upon the HST F702W image and 
the sub-mm conto ur map (solid thin black lines) presented in 
Borys et alJ (|2004 . Note that the SCUBA beam is 15 x 15 arc- 
sec, significantly larger than the VLA beam (6.99 x 6.03 arc- 
sec). To compare the radio and sub-mm emission at the same 
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Table 1. Details of the radio sources observed in the core of MS045 1.6-0305. The columns show: position (RA, DEC), peak 
flux density (Sp^^), total flux density (St) and deconvolved Gaussian sizes (major axis, minor axis and position angle) with their 
corresponding formal errors. Cases where the parameters of the Gaussian fits are not well constrained are indicated by a dash. 
The coordinates ar e given as off'sets with respect to the cluster centre, RA(J2000)=04:54:10.8 and DEC(J2000)=-03:00;51.6 (see 
Takata et al ]|2003l table 2). A version of this table in absolute coordinates can be found in the online material. 



Name 


RA 


DEC 




Sr 


Maj Axis 


Min Axis 


PA 




J2000 (") 


J2000 (") 




^Jy 


(") 


(") 


deg 


Ra2 


31.3±0.3 


1.1+0.4 


70 + 8 


95+ 18 


6+1 


2+1 


27 + 28 


Ral 


30.3±0.2 


-7.1+0.2 


109 + 9 


109 + 9 








Rbl 


26.5+0.2 


-26 + 0.2 


151+9 


151+9 








Rb2 


22.3±0.5 


-30.50 + 0.7 


52 + 8 


100 + 22 


9 + 2 


3 + 2 


158+14 


Rcl 


22.1±0.6 


-18.3 + 0.4 


50 + 9 


55 + 16 


6 + 2 




112+11 


Rc2 


13.5±0.6 


-31.3 + 0.8 


41 + 8 


58+18 


7 + 2 


1 + 3 


10+ 163 


Rc3 


35.4±0.8 


10.9 + 0.5 


52 + 8 


78+19 


8 + 2 


+ 2 


73 + 11 


Fa 


10.1+0.6 


-12.7 + 0.6 


45 + 9 


50+17 


4 + 6 


+ 4 


123 + 38 


Fb 


3.5 ±0.6 


-0.8 + 0.6 


49 + 9 


70 + 20 


6 + 5 


3 + 6 


122 + 45 


Fc 


50.4+0.6 


-37.5 + 0.6 


44 + 9 


44 + 9 








Fd 


64.1+0.1 


-17.2 + 0.1 


634 + 9 


1039 + 21 


7.6 + 0.1 


1.4 + 0.5 


128+1 


Fe 


-22.40+0.02 


-0.07 + 0.02 


1549 + 9 


1777 + 17 


3.3 + 0.1 


0.8 + 0.5 


125 + 3 



resolution (see Sect. |4]i, we also produced a tapered image of 
the radio data, weighting down the long baselines to reach a 
Gaussian restoring beam similar to the SCUBA beam. The re- 
sulting map is presented in Fig|3](solid white lines). 

The AIPS task IMFIT was used to fit Gaussian components 
to all the radio sources detected in the field (3 and 4 Gaussians 
simultaneously, in the case of the two extended regions of radio 
emission in the naturally weighted map). The radio positions 
obtained are represented by square boxes in FiglT] The results 
are listed in Table [T] together with their formal errors. 

Radio source components Ra, Rb, Rc and Fa appear to be 
related to the sub-mm emission. The source Fb may be related 
to the central brightest cluster galaxy (BCG) (see Table |6]l, but 
there is no obvious optical/NIR counterpart for Fc. 

We also detect two bright radio sources, Fd and Fe. The lat- 
ter is clearly identified with a optical/NIR counterpart, and both 
radio sources are almost two orders of magnitude brighter than 
the other radio sources in the field. These sources are proba- 
bly not lensed images of the same background source — Fe is 
more compact than Fd, even though it is brighter. The positions 
of Fd and Fe are coincide nt (within the errors) with two radio 
sources already reported in lStocke et al.l (119991) (see Table|6l). 



3. A lens model of the radio emission 

In order to aid our interpretation of the radio emission asso- 
ciated with SMM J04542-0301, we have created a simple el- 
liptical lens model with external shear to describe the cluster 
lens potential. Our an alysis employs the GRAVLENS software 
package developed by lKeetonl ( 2001 ). The modeling strategy is 
described in Sect l3.1l and the results of the model are discussed 
in Sectll2l 



3.1. Modeling strategy 

The first lens mo del of MS045 1 .6-0305 was presented in 
Takata et al.l (120031) . and describes the total mass distribution of 
the cluster by a singular isothermal ellipsoid. To be mor e sen- 
sitive to the local mass distribution, iBorvs et alj ( 120041) mod- 
eled the cluster core and 39 galaxy cluster members using 40 
smoothly truncated pseud o-isothermal ellipti cal mass distribu- 
tion profiles (PIEMD, see iKneib et al.lll996h . 

The critical curves of the best lens model found by 



Borvs et al.l (120041) are shown in FiglT] It can be seen that the 
tangential critical curve lies between the radio emission Ral- 
Ra2, Rbl-Rb2 and Rcl-Rc2. Based on the general properties of 
the lens geometry, this suggests that each of these image pairs 
belongs to a group of 3 images produced by one source located 
close to the caustic in the source plane. We propose the follow- 
ing scenario in an attempt to understand the radio emission we 
observe in terms of gravitational lensing: 

1 . Ra2-Ral are fold images of a source Ra, with an expected 
counterpart image {Ra3) close to ARC lei. 

2. Rbl-Rb2 are fold images of a source Rb with an expected 
counterpart image {Rb3) close to Rc3. 

3. Rcl-Rc2-Rc3 are multiple images of a single source Rc, lo- 
cated behind the cluster 

To test this hypothesis, we implemented a new lens model 
for MS045 1.6-0305 using the GRAVLENS code. Since we are 
only interested in testing the lensed nature of the radio emis- 
sion, we modeled the overall mass distribution of the cluster 
using a single mass profile plus external shear (see Appendix 
1 for a more detailed description of the related formulae). 
Unli ke the two previou s models, we choose an NEW pro- 
file dNavarro et al. I ll996h for the cluster mass distribution that 



is consistent with observations (see e.g. [ Pointecouteau et al 
20051: IComerford et alj|2006t iBassino et al.ll2006h and predic- 



tions from dark matter simulations. 
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To constrain the model, we performed a number of distinct 
steps. First, the parameters of the mass model were chosen in 
order to r eproduce the gener al shape of the critical lines deter- 
mined bv lBorvs et alj (120041) (see FiglTJ. Second, the positions 
and fluxes of the ERO images (Bl, CI, B2, C2, B3/C3) were 
used as constraints for the first optimisation of the input model. 
And third, we included the positions and fluxes of the radio 
images as new constraints (following the previous hypothesis) 
to re-optimise the model. The constraints are liste d in Tables [1] 
and|2 The coordinates of the cluster centre used in lTakata et al.l 
(l2003b were chosen as the origin of the coordinate system. 

The lens model obtained through this process turned out 
to have a degeneracy between the mass and the scale radius. 
To break this degeneracy, we used information about the con- 
centration parameter {6c) derived from ACDM N-body simu- 
lations. First, we produced a set of 10 new models varying the 
core radius between 40 and 150 (a range that contains the core 
radius value of the dege nerate model). Following the formalism 
presented in lBuUock et al. ( 2001 ). we calculated the concentra- 
tion parameter and the virial mass for the set of new models 
(see Appendix 1 for details). We found that, while 5c varies 
between 3.3 and 9.6, the virial mass is always of the order of 
IO'^Mq. For a hal o of that mass situat ed at z=0.55, the toy 
model presented in Bullock et al.l (120011) predicts a concentra- 
tion parameter of dc - 3.35. So we plotted 6c versus the core 
radius for the set of new models, interpolating the results in or- 
der to determine the core radius that corresponds to 6c = 3.35. 
Finally, we fixed the core radius of the model to this value and 
re-optimized the remaining parameters. 



Table 2. Optical/NIR constraints used in the lens model. The 

columns show: position (RA, DEC), total flux density in K' 
band and predicted magnification ju (see also Table 1 in B04). 
The coordinates are given as offsets with respect to the cluster 
cent re, RA(J2000)^04 :54:10.8 and DEC(J2000)=-03:00:51.6 
(see lTakata etal]|2003l Table 2). 



Name 


RA 


DEC 


Flux in K' band 








J2000 (") 


J2000 (") 








Bl 


29.95 


7.498 


3.6 ±0.1 


8± 


1 


B2 


28.47 


-9.702 


1.9 ±0.1 


10 ± 


1 


CI 


29.95 


4 


1.4 ± 0.1 


10 ± 


1 


C2 


28.47 


-7.502 


0.9 ±0.1 


5± 


1 


B3/C3 


1.481 


-32.9 


2.5 ±0.1 


5± 


1 



3.2. Results 

Figure |2] and Table |3] show the results of the lens model. 
Although the model is not unique, it is able to reproduce the 
positions of the ERO images and the radio emission (Ral, Ra2, 
Rbl, and Rb2) very well. The largest offsets are found for Rcl, 
Rc2 and Rc3 (see panel 4 in Fig|2]i which are the most distant 
images from the critical curves. This is consistent with the ef- 
fect of degeneracies in the global mass model near the critical 
curves: a change in the model parameters produces a more sig- 



nificant change in the image properties when they are located 
further from the critical curves. 

To improve the fit of images Rcl-Rc2-Rc3 , the redshift of 
the model was changed for source Rc, but this only produces 
a radial shift of their predicted positions iti the same direc- 
tion, something that cannot improve the fit shown in panel 4. 
Therefore, we believe that the most important contribution to 
this offset is probably coming from the group of galaxies in the 
region between Ral and Rcl (see FiglTJ. This group of galax- 
ies is expected to introduce perturbations in the overall mass 
distribution of the cluster which are not accounted for in the 
smooth NEW mass model. 

One notable result is that the model predicts two faint 
counterpart images {Ra3 and Rb3) that do not appear in our 
radio image. In our model Ra3 and Rb3 are less magnified 
than Ral-Ra2 and Rbl-Rb2 respectively. The predicted rela- 
tive magnifications suggest that Ra3 and Rb3 should appear in 
our maps at the 2 and 4 cr level respectively, but there is no ev- 
idence for this in the radio images. We note, however, that the 
predicted magnifications depend strongly on the overall mass 
model employed and perturbations by individual galaxies, so 
they should be treated as rough estimates of the true magni- 
fication. Therefore, the non-detection of these images do not 
necessarily mean that the lens model is wrong, since the real 
magnification could be less than that predicted by this simple 
model. 

The mass model parameters that characterise the NEW pro- 
file are summarised in Table |4] The errors represent the Icr 
level of the function of each parameter We note that the 
shear of the model is quite large, and the model centre is shifted 
ARA=-0.9"and ADEC=4.5" from the assumed position of the 
cluster centre. These effects are most likely systematic errors 
that compensate for the contribution of the group of galaxies 
that we are not including in the model, and the fact that we are 
forcing it to fit Rbl and Rb2 as mirror images (a hypothesis that 
we will discuss in detail in Sect|4]i. 

As already noted, this lens model is only meant to test if 
the configuration of the observed radio emission can be under- 
stood as the result of gravitational lensing. In that sense, Fig|2] 
shows that the model is able to reproduce the position of the 
EROs and the radio observations reasonably well, explaining 
the morphology of the radio map as the result of three lensed 
background radio sources. Therefore, we expect the two EROs 
and the three radio sources li e at the same redshi ft in the source 
plane (z ~ 2.9, estimated in iBorys et al.l (120041) for the EROs 



including the spectroscopic redshift information of the ARC in 
their lens model). 

But perhaps the most interesting result from the lens model 
is that the source Ra and the EROs B and C are predicted to be 
located inside a region of about 0.6 arcsec in the source plane 
(see panel 6 of Fig|2]i, which corresponds to a linear separation 
of only 4.7 kpc. The same situation is found in the case of the 
LBG and the ERO p air, which are sep arated by ~10 kpc in 
the source plane (see iBorys et al.ll2004l) . This means that the 
radio source Ra is lying just between the LBG and the EROs, 
all of which are located in a region smaller than the extent of 
a typical galaxy (~20 kpc). It would appear therefore, that the 
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Panel 1: EROs (Image Plane) Panel 2: Radio source Ra (Image Plane) 




Panel 3: Radio source Rb (Image Plane) 




# Rb2 



-40 -30 -20 -10 

arcsec 



Panel 5: Image Plane 



CD -10 





o Rc3 



-30 -20 



Panel 4: Radio source Rc (Imcge Plane) 

■ □Rc3 




Rc2 



-30 -20 



-10 







Of csec 
Panel 6: Source Plane 
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1^ Cluster's cen:er 
H BCG 
ERO B 

* ERO C 
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• Radio source Rb 
■ Radio source Rc 



-40 -30 -20 -10 

arcsec 



Fig. 2. Lens model results. The solid black curves represent the critical curv es (panels 1 - 5) an d caustics (panel 6) that define the 
model. Filled symbols represent the measured positions of the ERO images < Borvs et"al]|2Q04 ) and the radio components; empty 
symbols are the positions predicted by the model. The ERO pair and each group of suggested multiply-imaged radio components 
are shown seperately in different panels (1 to 4). Panel 5 shows all of these components plotted together. Panel 6 shows the 
position of the ERO pair and radio sources in the source plane. 
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Table 3. Lens model results. The columns show: coordinates of the predicted images (RA, DEC), predicted lensed flux (Sj-), 
offsets between the measured and predicted quantities (ARA, ADEC, ASj-) and predicted magnification (jj). The coordinates are 



given as off'sets with respect to the cluster centre, RA(J2000)=04:54:10.8 and DEC(J2000)=-03:00:51.6 (see iTakata et al .112003 
table 2). 



Name 


RA 


ARA 


DEC 


ADEC 


St 


A5t 






(") 


(") 


(") 


(") 


fj]y 


^Jy 




Bl 


-29.80 


-0.15 


7.25 


-2.75 


2.82 


0.75 


12.19 


B2 


-28.05 


-0.42 


-9.76 


0.06 


2.70 


0.83 


-11.66 


B3 


-1.50 


0.02 


-32.66 


-0.24 


1.13 


0.11 


4.89 


CI 


-30.07 


0.12 


4.23 


-0.23 


1.31 


0.11 


17.66 


C2 


-28.92 


0.45 


-7.40 


-0.1 


1.26 


0.36 


-17.00 


C3 


-1.47 


-0.01 


-33.14 


0.24 


0.35 


0.88 


4.75 


Ra2 


-31.03 


-0.29 


0.94 


0.17 


100.95 


-5.95 


25.18 


Ral 


-29.99 


-0.26 


-7.18 


0.12 


99.23 


-9.77 


-24.76 


Ra3 


-3.30 




-33.67 




19.73 




4.92 


Rbl 


-26.75 


0.23 


-26.52 


0.54 


131.12 


-19.88 


-36.78 


Rb2 


-22.11 


-0.22 


-30.70 


0.2 


111.05 


-11.05 


31.15 


Rb3 


-36.93 




-2.58 




37.65 




10.56 


Rcl 


-24.16 


2.11 


-22.12 


3.86 


72.17 


17.17 


-11.51 


Rc2 


-11.76 


-1.76 


-30.90 


-0.43 


62.34 


-4.34 


9.94 


Rc3 


-32.96 


-2.39 


10.28 


0.63 


42.37 


35.63 


6.75 



Table 4. Lens model parameters. The columns show: Mass scale (k^). Galaxy position in arcsec (xq, jq), ellipticity (e. Be), 
external shear (y, 0y) and scale radius in arcsec (r,). The errors represent the l-cr level of the fuction for each parameter 
The scale radius has no error estimations because its value was fixed during the optimization process to be consistent with a 
concentration parameter of 3.35 (see sect l3.1l l 



Mass Model 




yo e 




r 


e. 




NFW 0.331°-°\ 


-0.90!?;J 


4 51+0'' 55+""-'' 


-42.20!^-^ 


0.18t»^ 


61.44!^;:J 


[142.90] 



radio source Ra, the LBG and the ERO pair indeed constitute 
an interacting or merging system of galaxies. 

4. Comparison of the radio, sub-mm & optical/NIR 
data 

In this section we make a more detailed comparison between 
our r adio maps and the pre-existing sub-m m and NIR/optical 
data (iBorys et al.ll2004l:lTakata et al.l2003h . 



4.1. Sub-mm vs Radio emission: morphology and flux 
density ratio 

Figure |3] presents the sub-mm and tapered radio maps with a 
common resolution of ~ 15" x 15". The regions of radio and 
sub-mm emission are not only coincident, but they are extended 
on the same angular scale (~ 1') and have a strikingly similar 
morphology. This strongly suggests that the radio and sub-mm 
emission are associated with each other and are produced by 
the same (lensed) sources, probably star forming galaxies at z > 
2. In addition, the positions of the radio components located at 
Rbl, Rb2, Rcl, Rc2 and Fa (see FigUll are consistent with the 
sub-mm emission that coul d not be reproduce d using only the 
ERO images wAARCl (see lBorvs et"al]|2004 Fig.7). 



However, Fig|3] also shows some diiTerences in the mor- 
phology of the radio and sub-mm emission. The most relevant 
one is that the brightest region in sub-mm is not associated with 
the brightest region in radio {Rbl/Rb2), but with the second 
brightest {Ral/Ra2). One explanation for this apparent discrep- 
ancy in Rbl/Rb2, is that the radio and sub-mm emission arise 
from slightly different regions in the source plane and are dif- 
ferentially magnified. This effect could be quite significant for 
sources lying close to or extending across a caustic. Indeed, we 
note that recent Mid-IR and radio studies of local star form- 
ing galaxies show variations across the disk of u p to a factor of 



5 in the ratio of the FIR and radio luminosity dMurphv et al 



2004h . Another possibility is that the radio emission in the re- 



gion Rbl/Rb2 is not only associated with the sub-mm emission 
that arises from high-z star formation but from an additional 
component, perhaps an AGN in the foreground cluster that has 
no counterpart in the sub-mm. Indeed, a possible galaxy cluster 
member is located within 1 .5"of the radio component Rbl, and 
may be an optical/NIR counterpart to this source (see table|6]l. 

The tapered map (see Fig|3]l also shows an extension of the 
radio emission towards B3/C3 which is not seen in Fig[T] This 
suggests the possible existence of an extended radio source in 
this region, presumably associated with the faint su b-mm emis- 
sion "toe" that appears in the image presented by iBorvs et al 
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Fig. 3. The VLA 1 .36 GHz tapered contour map (solid white lines) superimposed upon the SCUB A 850-^m contour map (solid 
black lines) and the inverted HST F702W image of the centre of the cluster MS045 1 .6-0305 jBorvs et al.ll2004 . The axes 
represent the right ascension (x-axis) and declination (y-axis) in the J2000 coordinate system. Contours of the tapered radio map 
are drawn at -3, 3, 4, 5, 6, 8, 10, 11, 13 & 15 times the 1-cr noise level of 14.2 /^Jy/beam. Contours of the sub-mm map are drawn 
at4, 6,7, 9, 10, 11 & 1 1 .5 mJy/beam. The positions of the radio sources and relevant objects in the NIR (see fig. 1), are plotted as 
reference points. The black circle inside the box in the bottom-right corner is the beam size of the sub-mm map (15 x 15 arcsec), 
whereas the white one (almost covered by the black circle) corresponds to the beam-size of the radio map (15.06 x 14.26 arcsec, 
in position angle PA - 68.3°). 



(120041) . This is consistent with the existence of the radio coun- 
terpart Rc3 predicted by the lens model. 

Assuming that the radio and sub-mm emission is produced 
by the same galaxies, the 5'850/im/5' i.4ghz flux density ratio pro- 
vides information about their SEDs. Note that, since the sub- 
mm image does not resolve the radio components shown in 
FigHI due to the poor resolution of SCUBA, the radio flux 
density should be obtained from the tapered radio map. The 
integrated flux densities were calculated with the AIPS task 
TVSTAT, using the 4 mJy/beam sub-mm contour to delimit the 
same integration area in the sub-mm and tapered radio maps 
(see FigO. Using this method, we find^gso/im - 54.6 + 5.7mJy 
and S 1.4GHz = 0.547 ± 0.03mJy. The errors were calculated us- 
ing the expression (Tbeam x V^, where (Tbeam is the noise per 
beam of the image, and N is the number of beams within the 
area delimited by the 4 mJy/beam sub-mm contour. 

The observed S'gso^m/S' i.4ghz flux density ratio in 
MS' 045 1.6 - 0305 was compared with the flux density ratio 
obtained from a set of SED galaxy templates. This set is com- 
posed by the archetype star forming galaxies Arp220 and M82 
(Polletta et al. in prep), the AGN-dominated galaxy Mrk 231, 
an d the set of Blue Compact Dwarf galaxy SEDs presented 



m 



Hunt & Maiolinol (120051) . To detect possible diff'erences in 



the nature of faint and bright sub-mm sources, we also per- 
formed the same analysis for the faint sub-mm source detected 
in A2218, and the bright source detected in A1835 (see Table 
0. 



We find that the observed flux density ratio in MS 045 1.6- 
0305 is closer to the one obtained using the SED template of 
M82. However, we note that the observed total flux ratio is be- 
ing underestimated due to the "excess" of radio emission asso- 
ciated with Rbl/Rb2, and therefore a SED similar to Arp 220 
may be more appropriate for this source. Until the sub-mm 
emission can be resolved into diflFerent components to deter- 
mine their flux ratios independently, all we can conclude is that 
the overall flux density ratio is largely consistent with the fa- 
vored hypothesis that the bulk of the radio/sub-mm emission is 
arising from distant star forming galaxies that appear to follow 
the w ell known FIR-radio correlation ( Condon 1992: .Garrett! 
20021) . 

The bright source in A1835 is the only one which 
is well fitted by one o f the BCD templates presented by 
Hunt & Maiolinol (120051) . The faint source in A2218 seems 
to be similar to Mrk 231, perhaps suggesting it is a "warm" 
SCUBA source f ollowing the classification presented in 
Egami et al.l(l2004 . 



4.2. Optical/NIR counterparts to the radio emission 

Borvs et a D (|2004 proposed that the sub-mm emission is re- 
lated with three objects (see Fig[T]i: A LBG (imaged as ARCl 
and ARClci) and a pair of triply-imaged EROs {B, imaged as 
B1/B2/B3; and C, imaged as C1/C2/C3). 

As is shown in FiglT] ARCl is situated well inside the re- 
gion of radio emission associated with Ral/Ra2, and is there- 
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Table 5. Flux density ratios. The columns show: name of the sub-mm emission (SMM), redshift of the sub-mm source (z), 
cluster in which the sub-mm emision is located (Cluster lens), SED template that provides the best fit to the observed flux ratio 
(Best SED template), flux density ratio obtained from observations (FDR observed), flux densiy ratio predicted by the template 
(FDR template), deviation between the observed and predicted flux ratio (deviation). The value of FDR observed in SMM 
J16359+6612 was calculated by adding the fluxes of all the images together The deviations were calculated as (FDR template / 
FDR observed)-! . 



SMM z Cluster lens Best SED template FDR observed FDR template deviation References 



J1401 1+0252 


2.56 


A1835 


NGC5253 


127 + 


37 


115 


-0.09 


101 


J16359+6612 


2.516 


A2218 


Mrk231 


60 ± 


10 


49 


-0.18 


K04; G05 


J04542-0301 


2.9 


MS045 1.6-0305 


M82 


100 + 


11 


115 


0.15 


Sect. 4.1 



Table 6. Suggested NIR counterparts. The columns show: name of the radio source (Radio Source), name of the suggested 
NIR counterpart (Counterpart Source), coordinates of the counterpart source (RA cs, DEC cs), offsets between the position of 
the radio source and its NIR counterpart (ARA, ADEC) and references that contain information about the counterpart sources 
(References). Note that ARA, ADEC should be interpreted as indicative values (see Sect l4.2l for details). The coordinates are 



given as off'sets with respect to the cluster centre, RA(J2000)=04:54:10.8 and DEC(J2000)=-03:00:51.6 (see iTakata et al.ll2003 
Table 2). A version of this table in absolute coordinates can be found in the online material. 



Radio Source 


Counterpart Source 


RAcs 


DEC cs 


ARA 


ADEC 


References 






(") 


(") 


(") 


(") 




Fb 


BCG 


1.1 


-0.8 


2.4 


0.0 


Stockeetal. fl999) 


Fd 


0451-03C 


63.4 


-16.2 


0.8 


1 


Stockeetal. (1999) 


Fe 


045 1-03 A 


-21.6 


-1.2 


0.7 


1.1 


Stockeetal. (1999) 


Ral 


ARCl centre 


32.2 


-4.3 


1.9 


-2.8 


Borys, private comunication 


Ral 


ARCl bottom end 


31.9 


-6.3 


1.6 


-0.8 


Borvs et al. (2004). F720W HST image 


Ra2 


ARCl centre 


32.2 


-4.3 


0.8 


-3.2 


Borys, private comunication 


Ra2 


ARCl top end 


32.3 


-1.7 


0.9 


2.8 


Borvs et al. (2004). F720W HST imase 


Ra2 


Bl 


29.8 


-0.1 


-1.4 


-6.4 


Borvs et al. (2004) 


Ra2 


CI 


30.0 


0.1 


-1.4 


-2.9 


Borvs et al. (2004) 


Ral 


B2 


28.5 


-9.7 


-1.8 


2.1 


Borvs et al. (2004) 


Ral 


C2 


28.5 


-7.5 


-1.8 


0.4 


Borvs et al. (2004) 


Rbl 


Tc 


27.9 


-24.7 


1.4 


-1.3 


Takata et al. (2003) 


Rbl 


galaxy cluster member 


24.8 


-25.4 


-1.8 


-0.6 


Borvs et al. (2004), F720W HST image 


Rb2 


Td 


20.5 


-29.7 


-1.8 


-0.8 


Takata et al. (2003) 



fore probably related to it. However, in the case of the ERO 
pair, the images Bl and CI are located at the edge of the radio 
emission, suggesting that they are not directly contributing a 
significant amount of the radio flux density in this region. 

On the other han d, as shown in Fig [T] although ARClci 
and TF (an ERO from lTakata et al.ll2003b are most likely con- 
tributing to some of the measured flux density in the "sub-mm 
toe", B3/C3 is coincident with the maximum of this region. 
Therefore the EROs are certainly related to the sub-mm emis- 
sion, so we expect them to be related with the radio emission 
as well. 

The offsets between the estimated centres of the radio emis- 
sion and the optical/NIR candidates in the image plane are sum- 
marised in Tabled Note that the offsets of ARCl and the ERO 
images with respect to Ral and Ra2 are larger in declination 
(ADEC=3") than in right ascension (ARA=2"). This is proba- 
bly an effect of the magnification produced by the lens cluster, 
whose largest component is preferentially aligned in the direc- 
tion of declination (as is reflected in the direction of ARCl and 
the overall morphology of the sub-mm and radio emission). In 
the source plane, our lens model predicts that the offset between 



Ra and the ERO pair in the image plane is reduced to ~ 1". 
Indeed, as shown in Sect. 3. 2, the LBG and the EROs probably 
constitute an interacting system, with the radio source Ra situ- 
ated between them in the source plane. This suggests that the 
detected radio and sub-mm emission may come from the region 
in which the systems interact, perhaps due to the enhanced star 
formation produced by the merging process. A similar phenom- 
ena is also observed in the Antennae galaxy, where the bulk of 
the A2Q cm radio emission is situated between th e nuclei of 
both galaxies (see iHummel & van der Hulsl Il986l) . This sce- 
nario can explain the off'sets between the radio, sub-mm and 
NIR emission observed in the region of ARCl, Bl and CI in 
the image plane. 

Moving to the southern region of FiglU the open diamonds 
c orrespond to the pos itions of three additional EROs reported 



m iTakata et al.l (120031) . Two of these EROs {Tc and Td) are lo- 
cated within 2 arcseconds of Rbl and Rb2 (see table |6]l, and 
may be their NIR counterparts. In this case we expect them to 
be miiTor images, as we assumed in our lens model. We note 
that Takata et al. argue that these EROs have different photo- 
metric redshifts {zc - 3.730, zd - 0.5), which is inconsis- 
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tent with that hypothesis. However, both show the same colors 
(within the errors) in all bands except for B and 1^, and in those 
cases the differences may be due to contamination effects (from 
the galaxy cluster members situated close to Tc) and the use 



of a different aperture in each source (see iTakata et al.l 12003 



Table 1). Therefore, with this information we cannot discard 
the possibility that Tc and Td are lensed images of the same 
source. On the other hand, the scenario in which Tc and Td 
are not mirror images (but still the optical/NIR counterparts of 
Rbl and Rb2), is also possible within the lensing context. As is 
shown in Fi g[Tl Rbl is located very close to the critical curve 
presented in iBorys et al.l (120041) . so a small change in its po- 
sition can move it right on top of the critical curve (and the 
source component on top of the caustic), resulting in extremely 
high magnifications. This is consistent with the non-detection 
of the predicted counterpart image Rb3, and the high brightness 
of Rbl in the image plane. We also note that the shape of Tc is 
extremely elongated in the same direction of other faint arcs 
that appear in the same region of the Hubble image, suggesting 
that Tc may be lensed. 

Another possibility (as discussed earlier in SectlH, is that a 
possible cluster galaxy member could be the NIR counterpart 
of Rbl. This scenario can also explain the high brightness of 
Rbl if the galaxy cluster member turns out to be a "radio loud" 
AGN. 

Note that four of the six optical/NIR possible counterparts 
of the radio/sub-mm emission are EROs, wh i ch is c onsistent 
with the results presented in iKnudsen et al.l (120051) . On the 



other hand, studies carried out so far are incon clusive with re- 
spect to the overlap between LBGs and SMGs ( Chapman et al 



2002bHAdelberger et al.ll998HWebbl200ilHuang et al.l2005h 

However, in the scenario proposed here, the radio and sub-mm 
emission is the result of an interaction that involves a LBG, 
rather than emission coming from the LBG directly. 

Apart from the effects of lensing magnification, we also 
identify four other possible sources of error associated with the 
measured offset positions. In order of importance these include: 

- random measurement errors in the determination of the 
centre of the unresolved, blended radio components, 

- errors in the choice of the position of ARCl, due to its ex- 
tended and complicated structure, 

- systematic errors due to offsets between the HST and VLA 
coordinate reference systems (expected to be up to ~ 1", 
but clearly not dominant since no systematic trend is shown 
in Table |6]l. 

- intrinsic offsets that can appear if the radio and sub-mm 
emission come from different regions in the source plane. 

Since we cannot properly estimate the contribution of these 
errors (with the exception of the systematic error), the offsets 
shown in Table |6] should be taken as indicative values. 

5. Summary and Conclusions 

We have presented deep VLA archive observations at 1.4 
GHz of the central region of the cluster MS045 1.6-0305, 
discovering multiply-imaged radio counterparts to the sub- 
mm emission SMM J04542-0301, originally discovered by 



Chapman et al.l (I2002al) and recently studied by IBorys et al 
(l2004 . This is the second case of multiply-lensed radio 
emission coming from an intrinsically fai nt SMG (the first 
case was SMM J16 359-H6612 in A2218, see lKneib et alll2004 : 
Garrett et al.|[2005b . 



With a resolution of 7 x 6 arcseconds, the radio emission 
associated with SMM J0452-0301 can best be represented by 
seven discrete Gaussian components. A simple lens model 
of this system (based on a NFW mass profile) can reproduce 
the positions of the radio components assuming that they are 
multiple images of 3 background sources located at z-2.9. 
However, the model raises some questions that need to be 
resolved. 

Although the radio and sub-mm emission are clearly coin- 
cident and present a similar and unusually large angular extent 
(~1') and morphology (as expected if the radio and sub-mm 
emission comes from the same sources), the brightest peak of 
the radio emission is not coincident with the peak in the sub- 
mm. We find two possible scenarios that might explain this ob- 
servation: 

- the discrepancy is due to differential magnification pro- 
duced by the gravitational lensing effect - assuming that 
the radio and sub-mm emission have differnt morphologies 
and arise from different regions of the galaxy. Indeed we 
note, that R bl is situated very close to the critical curve 
presented in (iBorys et al ] l2004i see Fig.l). 

- the radio emission in that region includes a component as- 
sociated with an AGN associated with the foreground clus- 
ter. 

Borys et al. suggested that the sub-mm emission arises 
from an interacting system of galaxies formed by an ERO 
pair and a LBG. Although the association of the EROs with 
the radio emission appears to be uncertain in the image plane 
{Bl and CI are clearly located at the very edge of the radio 
emission), one of the three radio sources predicted by our lens 
model (Rd) is situated between the LBG and the ERO pair in 
the source plane. Our interpretation of this result is that the 
interacting region of the LBG and the ERO pair might be the 
source of the radio and sub-mm emission (due to the intense 
star formation generated during the merging process), whereas 
the optical/NIR emission might coiTespond to the cores of the 
merging galaxies. This scenario (a situation already observed 
e.g. in the Antennae galaxy) provides a consistent explanation 
of the offsets between the radio, sub-mm and NIR emission 
observed in the image plane (upper region of the map). 



^From the analysis presented in IBorys et alJ (120041) . it is 
also evident that the LBG and the ERO pair cannot account for 
all the emission coming from the central region of the sub-mm 
map. However, the higher resolution VLA observations show 
extended radio emission located in that region, which is 
expected to arise from 2 radio sources. Two of the components 
of that extended emission {Rbl and Rb2 ) seem to be related 
with another two EROs discovered by iTakata et alj ( l2003b 
{Tc and Td). Both EROs show similar colors, supporting 
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the idea (assumed in our lens model) that Rbl and Rb2 are 
images of the same source. However, this latter scenario 
im phes that the photometric redshifts of Tc and Td reported 
in iTakata et alj (|2003) may be incorrect. We also found a 
bright galaxy (probably a cluster member) that can be an 
optical/NlR counterpart of Rbl, keeping open the possibility 
that this particular region of the radio map might include 
emission from a foreground AGN. Unfortunately, none of the 
evidence is compelling enough to discriminate between the 
various scenarios that might explain the nature of Rbl and Rb2. 

In summary, we conclude that the radio and sub-mm 
emission found in MS045 1.6-0305 arises from at least 3 
highly magnified background so urces, one of them being the 
interacting system proposed by iBorvs et al.l (120041) (an LBG 
and an ERO pair). 

Further progress with this system requires a more complete 
comparison of multi-wavelength data to be made, and a more 
detailed lens model to be constructed. Deep, mid-IR observa- 
tions, as well as higher resolution sub-mm data, might be very 
important in understanding this system, in particular to confirm 
the possible lensed nature of Rbl and Rb2. We have recently 
re-observed MS045 1.6-0305 using the VLA at 1.4 GHz in its 
most extended A-array configuration. These higher resolution 
observations may shed new light on this system. 

Acknowledgements. The authors would like to thank Colin Borys for 
providing us with the HST/SCUBA reduced images and the con- 
tours of the lens model presented in his paper to prepare and analise 
our figures. We also want to thank James Bullock for providing 
us his code to estimate the concentration parameter expected for 
MS045 1.6-0305, Charles Keeton for helping to solve problems with 
the GRAVLENS code, Andy Biggs for his assistance during the data 
reduction, Antonio Heman Caballero and Leslie Hunt for providing 
us with the SED templates, and John Stockes for answering ques- 
tions related to the bright sources in the field. The authors are also 
very grateful to the referee for his/her constructive comments which 
helped to greatly improve the overall manuscript. ABA is also grate- 
ful to Edo Loenen for very valuable comments, suggestions and help 
during the writing of the manuscript. This work was supported by 
the European Community's Sixth Framework Marie Curie Research 
Training Network Programme, Contract No. MRTN-CT-2004-505183 
"ANGLES". 



References 

Adelberger, K. L., Steidel, C. C, GiavaHsco, M., et al. 1998, 
ApJ, 505, 18 

Bassino, L. R, Faifer, F R., Forte, J. C, et al. 2006, A&A, 451, 
789 

Blain, A. W. 1997, MNRAS, 290, 553 

Borys, C, Chapman, S., Donahue, M., et al. 2004, MNRAS, 

352, 759, (B04) 
Bryan, G. L. & Norman, M. L. 1998, ApJ, 495, 80 
Bullock, J. S., Kolatt, T. S., Sigad, Y., et al. 2001, MNRAS, 

321,559 

Chapman, S. C, Blain, A. W., Ivison, R. J., & Small, 1. R. 2003, 
Nature, 422, 695 



Chapman, S. C, Scott, D., Borys, C, & Fahlman, G. G. 2002a, 

MNRAS, 330, 92 
Chapman, S. C, Shapley, A., Steidel, C, & Windhorst, R. 

2002b, ApJ, 572, LI 
Comerford, J. M., Meneghetti, M., Bartelmann, M., & 

Schirmer, M. 2006, ApJ, 642, 39 
Condon, J. J. 1992, ARA&A, 30, 575 

Egami, E., Dole, H., Huang, J.-S., et al. 2004, ApJS, 154, 130 
Gaii-ett, M. A. 2002, A&A, 384, L19 

Gaii-ett, M. A., Knudsen, K. K., & van der Werf, R R 2005, 

A&A, 431,L21,(G05) 
Gioia, I. M. & Luppino, G. A. 1994, ApJS, 94, 583 
Huang, J.-S., Rigopoulou, D., Willner, S. R, et al. 2005, ApJ, 

634, 137 

Hummel, E. & van der Hulst, J. M. 1986, A&A, 155, 151 
Hunt, L. K. & MaioUno, R. 2005, ApJ, 626, L15 
Ivison, R. J., Greve, T R., Small, I., et al. 2002, MNRAS, 337, 
1 

Keeton, C. R. 2001, ArXiv Astrophysics e-prints, (astro- 
ph/0 102340) 

Kneib, J.-R, Ellis, R. S., Small, I., Couch, W. J., & Sharpies, 

R. M. 1996, ApJ, 471, 643 
Kneib, J.-R, van der Werf, R R, Kraiberg Knudsen, K., et al. 

2004, MNRAS, 349, 1211 
Knudsen, K. K. 2004, Ph.D. Thesis, leiden University 
Knudsen, K. K., van der Werf, R, Franx, M., et al. 2005, ApJ, 

632, L9 

Lilly, S. J., Bales, S. A., Gear, W. K. R, et al. 1999, ApJ, 518, 
641 

Luppino, G. A., Gioia, I. M., Hammer, F, Le Fevre, O., & 

Annis, J. A. 1999, A&AS, 136, 117, (L99) 
Murphy, E. J., Helou, G., Armus, L., et al. 2004, American 

Astronomical Society Meeting Abstracts, 205, 
Navarro, J. F, Frenk, C. S., & White, S. D. M. 1996, ApJ, 462, 

563 

Pointecouteau, E., Arnaud, M., & Pratt, G. W. 2005, A&A, 435, 
1 

Small, I., Ivison, R. J., & Blain, A. W. 1997, ApJ, 490, L5+ 
Smail, I., Ivison, R. J., Blain, A. W., & Kneib, J.-P 2002, 

MNRAS, 331,495 
Stocke, J. T, Perlman, E. S., Gioia, I. M., & Harvanek, M. 

1999, AJ, 117, 1967 
Swinbank, M., Small, I., Lindner, C, et al. 2006, ArXiv 

Astrophysics e-prints 
Takata, T., Kashikawa, N., Nakanishi, K., et al. 2003, PASJ, 55, 

789, (T03) 
Webb, T M. A. 2002, Ph.D. Thesis 



Appendix A: Calculation of the concentration 
parameter and the virial mass of the cluster 

The NFW density profile ( Navarro et alJ 19961) is defined in the 
three-dimensional space as: 



P(r) 



P.5 



(A.l) 



(r/r,)(l+r/r,)2 
where ps is the characteristic density and is the scale radius. 
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However the GRAVLENS code ( lKeetonll200ll) works with 
the projected surface mass density of the NFW profile, which 
is given by: 

\-T 

K(r) = Ik,— — - ; ^ PstJ^^va (A.2) 

X- - 1 

where x - rjr,, is the mass scale, and T is defined as: 



1 



= tan"' Vx2 - 1 (x > 1) 

■--\ 

' ^tanh"' Vl -x2 (jc < 1) 
(x= 1) 



Following the formalism used in iBuUock et alj (120011) . we 
define the scale radius as: 



(A.3) 



where c,,, is the concentration parameter and /?,,,> is the virial 
radius. 



Comparin g the definitions of virial mass u sed in 
Navarro et al.l d 19961) (M200) and iBullock et al] (I2OOII) (M,„>), 
the characteristic density can be written as: 



Ps - Puiz) 6c 



(A.4) 



where Pi,(z) is the universal density at redshift z, and Sc is the 
characteristic over-density, which is linked with c,,„ by the fol- 
lowing expression: 



f(Cvir) 



log(l+e,„v)--[|ii 



The parameter Ai, ,>(z) is called the virial o ver-density, and 
can be approximated dBrvan & Normanl[l998l) bv: 



A,„>(z) - (IStt^ + 82x - 39x2)/Q„,(z) - 1 
X = Q,„(z) - 1 



Combining equations lA.ll IA.2I IA.3I and I A. 41 the concen- 
tration parameter can be calculated using the following expres- 
sion: 

3 ^crit 



fiCrir) 



(A.5) 



p„(Zci) rs Ai,,>(Zci) 

The terms a-, and are given by the lens model (see ta- 
blegll. Using t he est imated redshift of MS045 1.6-0305 from 
Luppino et al. (Il999h (Zciustei—0-5 5) Pu(Zciuster) y^^s Ca lculated 



scaling the value for z=0 given in lBuUock et al.l ( 1200 lb (pu(z - 
0) = 8.3 X lO'^/i^MoMpc"^). 

Finally, knowing the value of c,,„ , the virial mass of the 
cluster can then be estimated as: 



47T 



-A,„>(zd)pu(zci) R; 



(A.6) 
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Table A.l. Details of the radio sources observed in the core of MS045 1.6-0305. The columns show: absolute coordinates in 
J2000 (RA,DEC), peak flux density (S/j^), total flux density (Sj-) and deconvolved Gaussian sizes (major axis, minor axis and 
position angle) with their corresponding formal errors. 



Name 


RA (+4'' 54"') 


DEC (-3°) 


s« 


Sr 


Maj Axis 


Min Axis 


PA 




J2000 (sec) 


J2000 (', ") 






ff 




deg 


Ra2 


12.89 ± 0.02 


00, 50.49 ± 0.43 


70 ± 8 


95± 18 


6± 1 


2± 1 


27 ±28 


Ral 


12.82 ± 0.01 


00, 58.66 ± 0.22 


109 ±9 


109 ±9 








Rbl 


12.57 ± 0.01 


01, 17.58 ±0.16 


151 ±9 


151 ±9 








Rb2 


12.29 ± 0.03 


01, 22.10 ±0.71 


52 ±8 


100 ± 22 


9±2 


3±2 


158 ± 14 


Rcl 


12.27 + 0.04 


01, 09.86 ±0.43 


50 ±9 


55± 16 


6±2 




112± 11 


Rc2 


1 1.70 ± 0.04 


01, 22.93 ±0.81 


41 ± 8 


58± 18 


7±2 


1 ±3 


10 ± 163 


Rc3 


13.16 ± 0.05 


00, 40.69 ± 0.46 


52 ±8 


78± 19 


8±2 


0±2 


73 ± 11 


Fa 


11.47 ±0.04 


01, 04.34 ±0.57 


45 ±9 


50 ± 17 


4±6 


0±4 


123 ±38 


Fb 


11.03 ±0.04 


00, 52.35 ± 0.59 


49 ±9 


70 ±20 


6±5 


3±6 


122 ± 45 


Fc 


14.16 ±0.04 


01, 29.06 ±0.55 


44 ±9 


44 ±9 








Fd 


15.08 ± 0.003 


01, 08.78 ±0.05 


634 ±9 


1039 ± 21 


7.61 ±0.13 


1.36 ±0.45 


128 ± 1 


Fe 


09.31 ±0.001 


00, 5 1.53 ±0.02 


1549 ± 9 


1777 ± 17 


3.3 ±0.1 


0.8 ±0.5 


125 ±3 



Table A.2. Suggested NIR counterpart sources of some radio emissions. The columns show: name of the radio source 
(Radio Source), name of the suggested NIR counterpart (Counterpart Source), absolute coordinates of the counterpart source 
in J2000 (RA cs, DEC cs), Radio-NIR offsets (ARA, ADEC) and references that contain information about the counterpart 
sources (References). 



Radio Source 


Counterpart Source 


RAcs (+4'' 54'") 


DECcs (-3°) 


ARA 


ADEC 


References 






(") 


(") 


(") 


(") 




Fb 


BCG 


10.73 


00, 50.81 


2.39 


0.04 


Stockeetal. (1999) 


Fd 


0451-03C 


15.04 


01,08.11 


0.75 


0.97 


Stockeetal. (1999) 


Fe 


045 1-03 A 


9.08 


00, 52.69 


0.73 


1.143 


^tocke et al. (1999) 


Ral 


ARCl centre 


12.946 


00, 55.900 


0.83 


1.9 


Borys, private comunication 


Ral 


ARCl bottom end 


12.926 


00, 57.900 


1.62 


-0.759 


Borvs et al. (2004)° 


Ra2 


ARCl centre 


12.946 


00, 55.900 


-3.19 


2.76 




Ra2 


ARCl top end 


12.953 


00, 53.280 


0.94 


2.79 


Borvs et al. (2004)° 


Rbl 


Tc 


12.66 


01, 16.3 


0.54 


1.28 


Takata et al. (2003) 


Rbl 


galaxy 


12.453 


01, 16.98 


-1.75 


-0.60 


Borvs et al. (2004)° 


Rb2 


Td 


12.17 


01, 21.3 


1.79 


0.80 


Takata, private communication 



from F720W HST image 



